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The  optimal  employment  of  autonomous  wide  area  search  munitions  is 
addressed.  The  scenario  considered  involves  an  airborne  munition  searching 
a  battle  space  for  stationary  targets  in  the  presence  of  false  targets.  Targets 
are  modelled  with  uniform,  Poisson,  and  normal  distributions.  False  targets 
are  modelled  with  Poisson  distributions.  All  relevant  parameters  can  be  ex¬ 
tracted  from  intelligence  information  on  the  enemy’s  order  of  battle  and  the 
sensor  performance  specification.  Analytic  weapon  effectiveness  measures 
are  derived  using  applied  probability  theory.  The  effectiveness  measures 
derived  in  this  paper  handle  time-varying  parameters  which  characterize 
the  battle  space  environment  and  the  performance  of  the  munition’s  sen¬ 
sor.  This  allows  the  formulation  and  solution  of  optimization  problems 
that  maximize  the  probability  of  a  target  attack  while  at  the  same  time 
constraining  the  probability  of  a  false  target  attack.  Optimal  schedules  for 
controlling  the  sensor  threshold  during  the  flight  are  derived  and  compared 
to  the  optimal  constant-threshold  results.  An  increase  in  weapon  effective¬ 
ness  is  demonstrated  when  the  sensor  threshold  is  dynamically  controlled 
during  the  flight. 
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I.  Introduction 


Several  types  of  autonomous  wide  area  search  munitions  are  currently  being  developed 
for  high  risk  air-to-ground  missions.  These  airborne  munitions  will  deploy  to  a  battle  space 
and  autonomously  search,  detect,  classify,  and  attack  targets.  In  this  paper,  the  presence  of 
false  targets  in  the  battle  space  is  acknowledged  np  front  and  analytic  weapon  effectiveness 
measures  are  derived  based  on  target  and  false  target  probability  distributions  as  well  as  the 
munition’s  sensor  characteristics.  All  relevant  parameters  can  be  extracted  from  intelligence 
information  on  the  enemy’s  order  of  battle  and  the  sensor  performance  specification.  This 
has  been  done  for  wide  area  search  munitions  in  the  past;  however,  the  problem  parameters 
were  assumed  constant.3  The  results  of  Ref.  3  are  now  generalized  to  handle  time-varying 
parameters.  This  then  allows  one  to  address  the  optimal  control  of  the  sensor  threshold  to 
maximize  the  probability  of  a  target  attack  without  unduly  increasing  the  probability  of  a 
false  target  attack.  The  latter  is  crucial  for  making  autonomous  operation  acceptable  to  the 
war  fighter. 

The  battle  space  geometries  considered  are  rectangular  and  circular.  In  a  rectangular 
battle  space  of  length  l  and  width  w,  the  munition  covers  the  area  using  a  constant  velocity 
v  starting  at  initial  time  t  —  0  and  ending  at  final  time  t  —  T.  The  sensor  footprint  is 
rectangular  with  swath  width  iv  and  incremental  swath  length  vdt  as  shown  in  Figure  1(a). 
Hence,  the  munition  covers  the  battle  space  in  one  sweep.  In  a  circular  battle  space  of 
radius  R,  the  munition  covers  the  area  using  concentric  annuli  of  thickness  dr  as  shown  in 
Figure  1(b).  The  munition  starts  at  the  origin  and  progresses  outward.  This  method  of 
concentric  annuli  approximates  an  outward  spiral  search  pattern. 


Figure  1.  Battle  Space  Geometries 
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The  scenarios  considered  involve  one  munition  searching  for  stationary  targets  among 
a  field  of  false  targets.  In  operations  research,  the  distribution  of  false  targets  in  a  battle 
space  is  often  modelled  using  a  Poisson  distribution.1  3  Three  such  scenarios  involving  a 
Poisson  field  of  false  targets  are  examined.  The  first  scenario  addresses  the  case  when  a 
single  target  may  appear  anywhere  within  a  rectangular  battle  space  with  equal  probability, 
thus  the  target  encounter  is  modelled  with  a  uniform  distribution.  The  second  scenario 
addresses  the  case  when  target  encounters  occur  at  an  expected  rate  in  a  rectangular  battle 
space,  thus  target  encounters  are  modelled  with  a  Poisson  distribution.  Finally,  the  third 
scenario  addresses  the  case  when  specific  information  is  available  about  a  target  (i.e. ,  position 
coordinates  with  some  error).  The  target  is  modelled  with  a  circular-normal  distribution  and 
the  corresponding  battle  space  is  defined  as  circular.  In  all  three  scenarios,  flight  altitude 
and  speed  are  constant,  and  the  munition  covers  the  entire  search  area. 

A.  Sensor  Performance 

When  a  sensor  encounters  an  object,  it  compares  the  image  to  a  stored  template  or  pattern 
and  either  declares  the  object  a  target  or  a  false  target.  In  practice,  detected  objects  are 
classified  to  a  certain  level  of  discrimination.  For  example,  one  may  classify  on  object  as 
either  air  breathing  or  ballistic.  A  finer  level  of  discrimination  may  be  a  specific  type  of 
air-breathing  or  ballistic  object.  Regardless  of  the  level  of  discrimination,  there  is  a  point 
where  the  sensor  reports  a  detected  object  as  either  a  target  thereby  authorizing  an  attack, 
or  a  false  target  thereby  commanding  no  attack.  Sensor  performance  is  judged  by  how  often 
the  sensor  is  correct.  The  probability  of  a  target  report,  Ptr,  is  the  probability  the  sensor 
correctly  reports  a  target  when  a  target  is  encountered.  The  probability  of  a  false  target 
report,  Pftr,  is  the  probability  the  sensor  correctly  reports  a  false  target  when  a  false  target 
is  encountered.  Together,  Ptr  and  Pftr  determine  the  entries  of  the  binary  “confusion 
matrix”  shown  in  Table  1  which  can  be  used  to  determine  the  outcome  of  a  random  draw 
each  time  an  object  is  encountered  in  simulation. 


Table  1.  Confusion  Matrix 


Declared  Object 

Encounl 

Target 

;ered  Object 
False  Target 

Target 

Ptr 

1  —  Pftr 

False  Target 

1  —  Ptr 

Pftr 

The  expression  (1  —  Ptr)  represents  the  probability  the  sensor  reports  a  false  target 
when  a  target  is  encountered.  This  type  of  error  results  in  a  target  not  being  attacked. 
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The  expression  (1  —  Pftr )  represents  the  probability  the  sensor  reports  a  target  when  a 
false  target  is  encountered.  This  type  of  error  results  in  a  false  target  being  attacked.  For 
this  binary  confusion  matrix,  true  positive  fraction  is  Ptr,  and  false  positive  fraction  is 
(1  —  Pftr)-  A  Receiver  Operating  Characteristic  (ROC)  curve  is  a  plot  of  true  positive 
fraction  versus  false  positive  fraction  that  starts  at  (0,0),  then  monotonically  increases  to 
(1, 1).  The  ROC  curve  model  adapted  from  Ref.  4  is  used: 


(1  —  Pftr) 


Ptr 

(1  -  c)  PTr  +  c 


(1) 


where  the  non-dimensional  scalar  parameter,  c  G  [l,oo),  depends  on  the  sensor  and  data 
processing  algorithm.  It  also  depends  on  the  munition  speed  (dwell  time),  and  engagement 
geometry,  which  includes  flight  altitude  and  look  angle.  A  family  of  ROC  curves  parame- 


Figure  2.  Family  of  ROC  Curves. 


terized  by  c  is  shown  in  Figure  2.  As  c  increases,  the  ROC  improves.  As  c  — >  oo,  the  area 
under  the  curve  approaches  unity  indicating  perfect  identification.  For  all  examples  in  this 
paper,  the  ROC  parameter  will  be  c  =  100. 

The  operating  point  on  a  ROC  curve  is  the  ordered  pair  [Ptr,  (1  —  Pftr)]-  For  a  given 
sensor  and  algorithm,  the  operating  point  is  determined  by  the  sensor’s  threshold.  Setting 
the  threshold  is  tantamount  to  setting  Ptr-  The  value  for  (1  —  Pftr)  is  calculated  using 
Eq.  (1).  Given  a  ROC  curve,  the  goal  is  finding  the  optimal  Ptr  -  either  constant,  or  as  a 
function  of  time  or  space. 
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B.  Optimal  Control  Problem  Formulation 

An  operator  wants  the  probability  of  a  target  attack,  PTa ,  in  a  battle  space  to  be  high  and 
the  probability  of  a  false  target  attack,  Pfta,  to  be  low.  Therefore,  Pta  and  Pfta  become 
weapon  effectiveness  measures.  Unfortunately,  increasing  Ptr  by  reducing  sensor  threshold 
with  a  view  to  increase  Pta  also  increases  (1  —  Pftr)  and  consequently  Pfta-  This  state 
of  affairs  is  modelled  by  the  ROC.  Thus,  the  operator’s  objectives  are  competing,  and  a 
trade-off  situation  arises.  To  ensure  Pfta  stays  low,  a  constraint  is  imposed  on  it.  The 
optimization  problem  statement  is  then 

Maximize:  Pta 
Subject  to:  Pfta  <  P ft Amax 

where  PFTAmax  is  set  by  the  designer  or  mission  planner. 

To  solve  the  posed  optimization  problem,  expressions  are  developed  for  Pta  and  Pfta 
for  the  considered  scenarios.  This  involves  integrating  probability  density  functions  parame¬ 
terized  by  Ptr,  velocity,  and  search  width,  as  well  as  data  on  the  battle  space  environment 
concerning  target  and  false  target  densities.  The  latter  are  obtained  from  knowledge  of  the 
order  of  battle.  Velocity  and  search  width  are  assumed  constant,  so  the  decision  variable 
is  Ptr-  The  explicit  ROC  model  given  by  Eq.  (1)  is  used;  however,  the  discerned  trends 
will  apply  to  any  ROC.  With  an  inequality  constraint,  the  unconstrained  problem  is  solved 
first  to  see  if  the  constraint  is  met.  If  Pfta  <  PFTAmax ,  die  unconstrained  solution  is  used. 
Otherwise  the  constraint  is  active  and  the  constrained  problem  is  solved  with  an  equality 
constraint.  Since  the  decision  variable  Ptr  is  a  probability,  it  is  bounded  by  [0, 1]. 

To  provide  a  benchmark,  constant-threshold  problems  are  solved  first.  Holding  Ptr 
constant  yields  a  static  (parameter)  optimization  problem  and  closed-form  solutions  are 
possible.  More  importantly  from  an  operational  point  of  view,  using  a  constant  sensor 
threshold  (i.e.  a  constant  Ptr)  is  the  easiest  way  to  employ  a  wide  area  search  munition.  The 
dynamic-threshold  problem  and  consequently  the  dynamic- Ptr  problem  can  be  formulated 
as  a  text  book  optimal  control  problem.5,6 

II.  Scenario  1:  One  Uniformly-Distributed  Target 
and  a  Poisson  Field  of  False  Targets 

The  scenario  when  one  target  is  uniformly  distributed  over  a  battle  space  among  a  Poisson 
held  of  false  targets  is  considered.  With  an  expendable  munition,  the  probability  of  an 
attack  occurring  during  the  time  interval  [t,  t  +  dt\  is  conditioned  on  the  probability  of  no 
attacks  occurring  prior  to  t.  The  probabilities  in  Ref.  3  were  calculated  assuming  constant 
parameters.  The  probabilities  are  re-derived  for  time-varying  parameters. 
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In  the  infinitesimal  interval  [t,t  +  dt],  like  that  shown  in  Figure  1(a),  three  events  can 
occur:  a  target  is  attacked,  a  false  target  is  attacked,  or  an  attack  does  not  occur.  Assuming 
no  attacks  occurred  prior  to  t,  the  probability  of  a  target  attack  occurring  in  the  time  interval 
[t,  t  +  dt]  is 

Ptr  (t)  ■  Pte  ( t )  (2) 

where  Pte  (t)  is  the  probability  of  a  target  encounter  in  the  infinitesimal  footprint.  Similarly, 
the  probability  of  a  false  target  attack  is 


[1  —  Pftr  (£)]  •  Pfte  (t)  (3) 

where  Pfte  (t)  is  the  probability  of  a  false  target  encounter  in  the  infinitesimal  footprint. 
Finally,  the  probability  of  no  attack  occurring  during  [t,  t  +  dt]  is 

[1  —  PTR  (t)]  •  Pte  (t)  +  Pftr  (t)  •  Pfte  (t)  •  (4) 


The  classification  probabilities,  Ptr  (t),  1  —  Ptr  (t),  Pftr  (t),  and  1  —  Pftr  ( t ),  are  specified. 
The  encounter  probabilities,  Pte  (t)  and  Pfte  (t),  are  determined  from  the  respective  target 
and  false  target  distributions.  For  one  uniformly-distributed  target, 


Pte  (t) 


w  (■ t )  v  (■ t )  dt 


A 


S 


(5) 


where  the  battle  space  or  area  searched,  As,  is 


A„  = 


w  (■ t )  v  (t)  dt. 


(6) 


For  a  Poisson  field  of  false  targets, 


Pfte  (t)  —  a(t)w  (t)  v  (t)  dt  (7) 

where  a  (t)  is  the  mean  false  target  density  at  time  t.  Indeed,  the  false  target  density  in  the 
battle  space  need  not  be  constant  and  could  vary  from  place  to  place  thereby  being  time 
dependent. 

Formulae  for  the  probabilities  of  no  previous  attacks  prior  to  t  are  now  given.  Detailed 
derivations  are  given  in  the  Appendix.  For  one  uniformly-distributed  target,  the  probability 
of  no  target  attack  occurring  prior  to  time  t  is 

pta  (t)  =  1  -  Ptr  (t)  w  (r)  v  (r)  dr.  (8) 
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For  a  Poisson  field  of  false  targets  with  density  a  (r),  0  <  r  <  t,  the  probability  of  no  false 
target  attacks  occurring  prior  to  time  t  is 


p  (t)  =  e~  /o[1_PFTil(r)]a(r)«)(rMT)(iT 


(9) 


The  probability  of  a  target  attack  can  now  be  calculated  as 


Pta  —  /  Ptr  (t)  ■  Pte  (t)  ■  Pfta  (t) 


—  Ptr  (t)  w  (t)  v  (; t )  e~  (10) 

'0 

and  the  probability  of  a  false  target  is 

Pfta  =  [  [1  —  Pftr  (t)]  •  Pfte  (t)  ■  Pta  (t)  '  Pfta  (t) 

Jo 

=  f  [1  -  Pftr  (t)]  a  (t)  w  (t)  v  (t)  x 
Jo 

^1  -  J  ~XPTR  (t)  W  (r)  V  (r)  dr\  e~  /o  [l-PFTR(T)}a{r)W{r)v{r)dr  dt^ 

The  problem  of  finding  the  optimal  Ptr  (t),  0  <  t  <  T,  assuming  w,  v,  and  a  are  constant 
is  now  addressed.  Recall  [1  —  Pftr  (£)]  is  determined  from  PTr  (t),  because  they  are  related 
by  the  ROC. 


A.  Constant-Threshold 


When  the  sensor  threshold  is  constant  during  the  flight,  Ptr  and  consequently  Pftr  are 
constant,  and  closed-form  solutions  are  obtained  when  evaluating  the  integrals  in  Eq.  (10) 
and  Eq.  (11).  The  probability  of  a  target  attack  is 


Pta 


Ptr 


l  _  g-(l -Pftr)oiwvT 

(1  —  Pftr )  cnwvT 


(12) 


and  the  probability  of  a  false  target  attack  is 


Pfta  -  >  -  'ta  -  M  -  Ptr)  e  "  . . 


(13) 


The  ROC  model  in  Eq.  (1)  is  invoked  to  reduce  the  problem  to  a  one-dimensional  search  for 
Ptr-  The  goal  is  maximizing  PTa  subject  to  Pfta  <  P FTAmax  ■ 

To  show  an  example,  values  are  needed  for  the  constants  c,  w,  v,  a,  and  T.  Units  are  time 
units.  Assuming  c  =  100,  aw v  =  50  [1/time]  and  T  =  0.5  [time],  the  outcome  probabilities 
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are  plotted  versus  Ptr-  Figure  3  shows  the  best  unconstrained  solution  is  P^n  =  0.723 
with  a  corresponding  PRA  =  0.535  and  Pftau  —  0.318.  The  subscript  u  is  used  to  denote 
unconstrained  solutions.  If  Pfta  is  bounded  by  PFTAmax  =  0.2,  the  best  constrained  solution 
is  PRR  =  0.563  with  a  corresponding  PRA  =  0.483.  In  Figure  3,  the  outcome  probability 
functions  are  smooth  and  well-behaved.  The  function  for  Pta  has  only  one  peak  and  never 
crosses  the  line  where  Pta  =  Ptr-  The  function  for  Pfta  is  monotonically  increasing  in 
Ptr,  so  any  constrained  solution  will  be  unique.  Finally,  the  function  for  the  probability  of 
no  attack  occurring  is  monotonically  decreasing  in  Ptr  as  expected. 


Figure  3.  Outcome  probabilities  versus  probability  of  a  target  report.  The  scenario  involves 
one  uniformly-distributed  target  among  a  Poisson  field  of  false  targets.  Probability  of  a  target 
report  is  constant  throughout  the  flight,  c  =  100,  awv  =  50  [1/time],  and  T  =  0.5  [time]. 


The  constant-threshold  optimization  problem  can  be  solved  for  a  number  of  PFTAmax  val¬ 
ues  and  a  plot  of  PRR  versus  P FTAmax  can  be  generated.  Figure  4  illustrates  the  sensitivity  of 
the  solution  to  changes  in  the  upper  bound  P FTAmax •  Tor  P FTAmax  >  0.32,  the  unconstrained 
solution  is  used. 
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Maximum  Allowable  Probability  of  a  False  Target  Attack 

Figure  4.  Optimal  probability  of  a  target  report  versus  maximum  allowable  probability  of  a 
false  target  attack.  The  scenario  involves  one  uniformly-distributed  target  among  a  Poisson 
field  of  false  targets.  Probability  of  a  target  report  is  constant  throughout  the  flight,  c  =  100, 
awv  =  50  [1/time],  and  T  =  0.5  [time]. 


B.  Dynamic-Threshold 

When  Ptr,  and  consequently  Pftr,  are  dynamic,  an  optimal  control  problem  is  on  hand. 
The  state  variables  become 


x  —  awv  [1  —  Pftr  (t)]  dr 


y=  -Ptr  (t)  dr 


z  = 


awv  [1  -  PFTr  (t)]  [1  -  y  (r)]  e  x(T)dr 


(14) 

(15) 

(16) 


and  the  cost  functional  integrand,  L ,  is 


L=--PTR(t)e-‘"K 


(17) 


Obviously, 


Pta  =  ~  L  [PTR  (t) ,  x  (t)]  dt. 


(18) 
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The  states  equations  are 


x  =  awv  [1  —  Pftr  (£)]  (19) 

V  =  t^Ptr  ( t )  (20) 

z  =  awv  [1  -  PFTr  (£)]  [1  -y  (£)]  (21) 

After  applying  Eq.  (1)  for  the  ROC  model,  the  Hamiltonian,  H ,  becomes 

H  =  (A,  -  e-)  1 PTR  +  awv— - [A,  +  A,(l  -  v)  (22) 

J  0  —  CJ  fTR  +  C 

and  the  costate  differential  equations  become 

Ax  =  A zawv— - (1  -  y)  e~x  -  ^PTRe~x  (23) 

(1  -  c)  PTR  +  c  T 

A  y  =  XzawvT - — — e-*  (24) 

(1  -  c)  PTR  +  c 

A,  =  0  (25) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  PRR  gives 


dH  _  ,  _x\  1  [Xx  +  A2  (1  -  y)  e  x]  awvc 

dPTR~[v  6  ]T+  [(1  -C)PTR  +  Cf 


Solving  =  0,  the  optimal  control  is 


c  ± 


Ptr  (t)  — 


{Ax(t)+A*[l-j/(t)]e  *W}atui; cT 


e  x(t)  _ \y(t) 


(27) 


Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  PRR  outside  of  [0, 1]  for  all 
time.  One  is  also  interested  in  the  optimal  unconstrained  solution,  PRRu ,  which  is 


Ptru  (t) 


c  —  \J Xx  ( t )  awvcTex^> 
c  —  1 


(28) 


Taking  the  derivative  of  Eq.  (28)  with  respect  to  t  gives 


PtRu  (t)  ~ 


1  aivve 


-x(t) 


4cTAx  ( t ) 


[Ptr  (0]  >  0- 


(29) 
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Therefore,  Prru  is  monotonically  increasing.  At  the  terminal  point  where  Xx  (T)  =  0, 
Ptru  (T)  =  ~Zi  >  1-  Thus,  the  constraint  Ptr  G  [0, 1]  becomes  active  prior  to  the  end 
and  remains  active  until  time  T.  This  represents,  for  the  unconstrained  problem,  a  “go  for 
broke”  tactic  in  the  end  game. 

Assuming  c  =  100,  awv  =  50  [1/time],  and  T  =  0.5  [time],  the  constrained  problem  can 
be  solved  for  a  number  of  P FTAmax  values.  Figure  5  shows  optimal  solutions  as  a  function 
of  time  for  three  values  of  PFTAmax  •  As  P FTAmax  decreases,  the  optimal  solution  flattens  out 
approaching  a  constant-threshold  solution. 


Figure  5.  Optimal  probability  of  a  target  report  versus  time.  The  scenario  involves  one 
uniformly-distributed  target  among  a  Poisson  field  of  false  targets.  c=  100,  awv  =  50  [1/time], 
and  T  =  0.5  [time]. 


The  objective  function  values  of  the  dynamic  solutions  can  be  compared  to  those  of  the 
constant  solutions.  Figure  6  shows,  for  the  assumed  parameters,  that  the  optimal  dynamic 
solution  is  at  most  3.4  percent  better  than  the  optimal  constant  solution.  Depending  on  the 
application,  this  improvement  may  or  may  not  be  worth  the  added  complexity. 
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Maximum  Allowable  Probability  of  a  False  Target  Attack 

Figure  6.  Optimal  probability  of  a  target  attack  versus  maximum  allowable  probability  of  a 
false  target  attack.  The  scenario  involves  one  uniformly-distributed  target  among  a  Poisson 
field  of  false  targets,  c  =  100,  awv  =  50  [1/time],  and  T  =  0.5  [time]. 


III.  Scenario  2:  Poisson  Fields  of  Targets  and  False  Targets 

The  scenario  when  targets  are  also  modelled  with  a  Poisson  distribution  is  considered. 
For  this  scenario  the  battle  space  consists  of  a  Poisson  field  of  targets  and  a  Poisson  field 
of  false  targets.  The  derivation  of  probabilities  is  similar  to  the  previous  scenario  with  the 
following  changes.  For  a  Poisson  field  of  targets, 

Pte  (t)  —  (3  (t)w  (t)  v  (t)  dt  (30) 

where  f3  (t)  is  the  mean  target  density  at  time  t.  The  probability  of  no  target  attacks 
occurring  prior  to  time  t  is 


p—  =  g-  fo  )w(t )v(r)dr 


(31) 


The  probability  of  a  target  attack  can  now  be  calculated  as 


Pta  —  /  Ptr  (t)  *  Pte  (t)  •  Pta  W  *  Pfta  W 


Ptr  (' t )  f3  ( t )  w  (; t )  v  (: t )  e~  ^ 
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and  the  probability  of  a  false  target  attack  is 


P 


FTA 


[1  -  PfTR  (£)]  •  PfTE  (t)  ■  Pta  (t)  ■  P~FTA  (t) 


[1  -  PfTR  (t)]  a  (: t )  W  (t)  V  (t)  e-Io{^-pFTR(T)nr)+PTR(T)0(r)}W(rMr)drdt  (33) 


The  parameters  w,  v,  a,  and  (3  are  assumed  constant. 


A.  Constant-Threshold 


When  Ptr  and  consequently  Pftr  are  constant,  closed-form  solutions  are  obtained  when 
evaluating  the  integrals  in  Eq.  (32)  and  Eq.  (33).  The  probability  of  a  target  attack  is 

p  —  _ PTR  ft _  _  e-[(l-PFTR)ot+PTR&}wvT 

(1  —  Pftr )  a  +  Ptr(3 
and  the  probability  of  a  false  target  attack  is 


-  1  pm  .  _  ^-[(1  -Pftr)oi+Ptr(3}wvT 

rFT A  ~  1  TfA  e 


(35) 


To  show  an  example,  values  must  be  provided  for  the  various  constants.  Units  are  time  and 
distance  units.  Invoking  the  ROC  model  in  Eq.  (1)  and  assuming  c  =  100,  w  =  0.2  [distance], 
v  =  50  [distance/time],  a  =  10  [1/distance2],  (3  =  1  [1/distance2],  and  T  =  0.5  [time], 
the  outcome  probabilities  are  plotted  versus  Ptr-  Figure  7  shows  the  best  unconstrained 
solution  is  Ptru  —  0.513  with  a  corresponding  Pta„  =  0.793  and  Pftau  =  0.161.  If  Pfta  is 
bound  by  P FTAmax  =  0.1,  the  best  constrained  solution  is  P^R  =  0.292  with  a  corresponding 
P^A  =  0.711. 

Like  the  previous  scenario,  the  function  for  Pta  has  only  one  peak.  Unlike  the  previous 
scenario,  the  values  of  a  and  (3  can  skew  the  function  for  Pta  above  or  below  the  line  where 
Pta  =  Ptr-  However,  the  function  for  Pfta  is  still  monotonically  increasing.  Thus,  any 
constrained  solution  will  be  unique. 


B.  Dynamic-Threshold 

When  Ptr,  and  consequently  Pftr,  are  dynamic,  the  state  variables  become 

x  =  f  {[  1  —  P ftr  (t)]  a  +  Ptr  (t)  (3}  wvdr  (36) 

Jo 

z  =  f  [1  —  Pftr  (t)]  awve~x^ dr  (37) 

Jo 
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Figure  7.  Outcome  probabilities  versus  probability  of  a  target  report.  The  scenario  involves 
a  Poisson  field  of  targets  and  a  Poisson  held  of  false  targets.  Probability  of  a  target  report 
is  constant  throughout  the  flight,  c  =  100,  w  =  0.2  [distance],  v  =  50  [distance/time],  a  =  10 
[1/distance2],  (3  =  1  [1 /distance2],  and  T  =  0.5  [time]. 


and  L  is 

L  =  -PTR(t)  (3wve~x{t).  (38) 

Since  both  targets  and  false  targets  have  Poisson  distributions,  the  state  x  represents  a 
combined  dynamic  Poisson  parameter.  Hence  there  is  no  need  for  a  separate  y  state.  The 
state  equations  are 


x  =  { [1  -  PFTr  (£)]  a  +  PTr  (t)  /?}  wv 
z  =  [1  —  Pftr  (£)]  awve~x® 

After  applying  Eq.  (1)  for  the  ROC  model,  the  Hamiltonian  becomes 

H  =  (Xx-  e~x)  PTr/3wv  +  (Ax  +  Xze~x )  awv- - ^77 — — 

(1  -c)  PTr  +  c 

and  the  costate  differential  equations  become 


A,  = 


A. .ft 


P, 


tr 


(1  -  c)  Ptr  +  c 


—  PtrP 


zuve 


(39) 

(40) 


(41) 


(42) 

(43) 
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Therefore,  Ptr,,  is  monotonically  increasing  and  PRR  (T)  =  Ay  >  1.  Thus,  the  “go  for 
broke”  tactic  appears  again  for  this  unconstrained  problem. 

Assuming  c  =  100,  w  =  0.2  [distance],  v  =  50  [distance/time],  a  =  10  [1/distance2] , 
(3  =  1  [1/distance2],  and  T  =  0.5  [time],  the  constrained  problem  is  solved  for  a  number  of 
PFTAmax  values.  Figure  8  shows  optimal  solutions  as  a  function  of  time  for  three  values  of 
P FTAmax  ■  Once  again,  as  P FTAmax  decreases,  the  optimal  solution  flattens  out  approaching  a 
constant-threshold  solution. 

The  objective  function  values  of  the  dynamic  solutions  are  compared  to  those  of  the 
constant  solutions.  Figure  9  shows,  for  the  assumed  parameters,  that  the  optimal  dynamic 
solution  is  at  most  3.6  percent  better  than  the  optimal  constant  solution.  Depending  on 
the  application,  this  improvement  may  or  may  not  be  worth  the  added  complexity.  These 
first  two  scenarios  involved  constant  ratios  of  targets  to  false  targets,  regardless  of  the  sensor 
footprint  position.  The  next  scenario  involves  non-constant  ratios  of  targets  to  false  targets. 
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Figure  8.  Optimal  probability  of  a  target  report  versus  time.  The  scenario  involves  a  Pois¬ 
son  field  of  targets  and  a  Poisson  field  of  false  targets,  c  =  100,  w  =  0.2  [distance],  v  =  50 
[distance/time],  a  =  10  [1/distance2],  (3  =  1  [1/distance2],  and  T  =  0.5  [time]. 


Maximum  Allowable  Probability  of  a  False  Target  Attack 

Figure  9.  Optimal  probability  of  a  target  attack  versus  maximum  allowable  probability  of 
a  false  target  attack.  The  scenario  involves  a  Poisson  field  of  targets  and  a  Poisson  field  of 
false  targets,  c  =  100,  w  =  0.2  [distance],  v  =  50  [distance/time],  a  =  10  [1/distance2],  (3  =  1 
[1 /distance2],  and  T  =  0.5  [time]. 
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IV.  Scenario  3:  One  Normally-Distributed  Target 
and  a  Poisson  Field  of  False  Targets 


The  scenario  when  one  target  is  normally  distributed  over  a  circular  battle  space  among 
a  Poisson  field  of  false  targets  is  now  considered.  The  ratio  of  targets  to  false  targets  depends 
on  the  sensor  footprint  position.  The  normal  distribution  is  circular,  centered  at  the  origin 
with  variance  a2.  The  derivation  of  probabilities  is  similar  to  the  first  scenario  with  the 
following  changes.  For  one  normally-distributed  target, 

1  _  r2 

Pte  (r)  —  - — 2^27r rdr.  (48) 

2no- 

The  probability  of  no  false  target  attacks  occurring  prior  to  reaching  radius  r  is 


Pfta  (r) 


—  e~  fo  2-n[l-pFTR(p)]a(p)pdp 


(49) 


and  the  probability  of  no  target  attack  occurring  prior  to  reaching  radius  r  is 


Pta  (r)  —  1 


Ptr(p)  ^ 
- 2 — e  2^  pdp. 


Jo  & 

The  probability  of  a  target  attack  can  now  be  calculated  as 

fR 

Pta  —  /  Ptr  (t)  ■  PtfAt)  • 


fR. 


Ptr  (r)  —,e-^  •  e"  J'o^H-PFTihP)] a{P)PdPdr 
a2 


and  the  probability  of  a  false  target  attack  is 


fR 


Pfta  —  /  [1  —  Pftr  (?")]  •  Pfte  (r)  ■  Ppj  (r)  ■  Pfta  (r) 


fR 


[1  -  PFTR  (r)]  27rra  (r)  x 


1  -  /  Ptr  (p)  -^e  ^  dp 

./  n  <7 


g-  fo  2tt[1  -PFTR(p)]a(p)pdp 


(50) 


(51) 


(52) 


The  parameter  a  is  assumed  constant. 
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A.  Constant-Threshold 


If  Ptr  and  consequently  Pftr  are  constant,  closed-form  solutions  are  obtained  when  evalu¬ 
ating  the  integrals  in  Eq.  (51)  and  Eq.  (52).  The  probability  of  a  target  attack  is 


Pta  — 


P 


TR 


2a2  [(1  -  Pftr)  om  +  ^ 


^1  -  q-^-Pftr^+^r2  j 


and  the  probability  of  a  false  target  attack  is 


(53) 


P ft  a  =  1  —  Pta  ~  1  —  P; 


tr 


_ 

+  Ptrc^2 


3-(l -PFtr)qfR2 


(54) 


To  show  an  example,  values  must  be  provided  for  the  various  constants.  Units  are  distance 
units.  Invoking  the  ROC  model  in  Eq.  (1)  and  assuming  c  =  100,  a  =  5  [1/distance2], 
a2  =  0.25  [distance2],  and  R  —  5  [distance],  the  outcome  probabilities  are  plotted  versus 
Ptr ■  Figure  10  shows  the  best  unconstrained  solution  is  PrR  =  0.788  with  a  corresponding 
Pta.u  =  0-615  and  Prtau  =  0.385.  If  Pfta  is  bound  by  P FTAmax  =  0.20,  the  best  constrained 
solution  is  PRR  =  0.057  with  a  corresponding  PRA  =  0.057.  This  is  a  significant  reduction  in 
performance. 


Figure  10.  Outcome  probabilities  versus  probability  of  a  target  report.  The  scenario  involves 
a  normally-distributed  target  and  a  Poisson  field  of  false  targets.  Probability  of  a  target  report 
is  constant  throughout  the  flight,  c  =  100,  a  =  5  [1/distance2],  a2  =  0.25  [distance2],  and  R  =  5 
[distance] . 


Like  the  scenario  with  one  uniformly-distributed  target,  the  function  for  Pta  has  one 
peak  and  never  crosses  the  line  where  Pta  =  Ptr ■  However,  as  a  result  of  the  changing 
ratio  of  targets  to  false  targets,  the  function  for  Pfta  is  not  necessarily  increasing.  Figure  10 
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shows  Pfta  increases,  then  decreases,  then  increases  again.  For  values  of  Pfta  above  the 
local  minimum  at  Pfta  =  0.385,  the  unconstrained  solution  is  used.  On  a  plot  of  PRR  versus 
P FTAmax  j  a  discontinuous  jump  appears  as  shown  in  Figure  11. 


Maximum  Allowable  Probability  of  a  False  Target  Attack 

Figure  11.  Optimal  probability  of  a  target  report  versus  maximum  allowable  probability  of  a 
false  target  attack.  The  scenario  involves  a  normally-distributed  target  and  a  Poisson  field  of 
false  targets.  Probability  of  a  target  report  is  constant  throughout  the  flight,  c  =  100,  a  =  5 
[1/distance2],  a2  =  0.25  [distance2],  and  R  =  5  [distance]. 


B.  Dynamic-Threshold 

When  Ptr,  and  consequently  Pftr ,  are  dynamic,  the  state  variables  become 


x 

y 

z 


and  L  is 


The  state  equations  are 


27 rap  [1  -  Pftr  (p)]  dp 

Ptr  (p)  — re  dp 
a 2 

2nap  [1  -  Pftr  (p)]  [1  ~V  (p)]  e~x{p)dp 


L  —  —Prr  (r) 


—e-^e~x{r) 


x  =  27 rar  [1  -  Pftr  (r)] 

r  _  r2 

y  =  Ptr  ( r ) 

a 2 

=  2 imr  [1  -  Pftr  (r)]  [1  -  y  (r)]  e-'1’'1 


(55) 

(56) 

(57) 


(58) 


z 
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(59) 

(60) 
(61) 


After  applying  Eq.  (1)  for  the  ROC  model,  the  Hamiltonian  becomes 


H=  (Xy-e  x)  PTr  (r)  ^  e  +  2vr ar 


TR 


(1  -  c)  PTR  +  C 


Ax  +  Xz  (1  —  y)  e  x 


and  the  costate  differential  equations  become 

•  Ptr 

Ac  =  A227TW— - — - 

(1  -  c)  PTR  +  C 

Ptr 

A  y  =  \z2nar— - — - —  t 

(1  -  c)  PTR  +  c 

A,  =  0 


(1  ~y)e  x  -  PTR(r)  —e  e 


(62) 


(63) 

(64) 

(65) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 

9H  _(  x  r  -A  |  [K  +  Xz(l-y)e-x]2narc 

oPtr  <7  [(1  —  c)  Ptr  +  cj 

Solving  =  0,  the  optimal  control  is 


c  ± 


Ptr  (r) 


1  [Aa;  (r)+Az  [1 

,  r2 

— y  (r )]  e ~  x  (r )  J  27raca2  e  2°^ 

/ 

e  x(r)—\ y(r) 

c-  1 


(66) 


(67) 


Once  again,  only  the  “minus”  root  is  used.  One  is  also  interested  in  the  optimal  unconstrained 
solution,  which  is 


p* 

rTRu 


r  = 


Xx  (r)  27raca2e 


*(»■)+£' 


Taking  the  derivative  of  Eq.  (68)  with  respect  to  r  gives 


(68) 


P 


TRU 


r  = 


2nacr2 


\  4(c-l); 


a2  At 


e[*W+^' 


Ax  (r) 


i) 


•  (69) 


Unlike  the  previous  scenarios,  it  is  unclear  whether  Ptr,,  is  increasing  or  decreasing. 

Assuming  c  =  100,  a  =  5  [1/distance2],  a2  =  0.25  [distance2],  and  R  =  5  [distance],  the 
constrained  problem  is  solved  for  a  number  of  P FTAmax  values.  Figure  12  shows  as  P FTAmax 
decreases,  the  sensor  is  “turned  off’  earlier  in  the  search. 
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Figure  12.  Optimal  probability  of  a  target  report  versus  radius.  The  scenario  involves  a 
normally-distributed  target  and  a  Poisson  field  of  false  targets,  c  =  100,  a  =  5  [1/distance2], 
a2  =  0.25  [distance2],  and  R  =  5  [distance]. 


The  objective  function  values  of  the  dynamic  solutions  are  compared  to  those  of  the 
constant  solutions.  Figure  13  shows  the  poor  performance  of  the  constant-threshold  solution 
when  PFTAmaX  <  0.385.  Unlike  the  previous  two  scenarios  where  dynamic-threshold  solu¬ 
tions  were  only  a  few  percent  better  than  constant-threshold  solutions,  there  is  substantial 
improvements  when  PFTAmax  <  0.385.  Since  the  dynamic  problem  allows  the  sensor  to  be 
“turned  off’  at  some  point,  the  solutions  are  similar  to  those  at  smaller  R  values.  In  essence, 
one  learns  how  far  out  in  radius  to  search.  The  dynamic  problem  “recognizes”  the  benefit 
of  staying  close  to  the  origin  when  a  target  is  normally  distributed.  The  constant-parameter 
problem  forces  one  to  arbitrarily  pick  R  then  optimize  Ptr  based  on  R.  Thus,  one  must 
settle  for  poor  performance  when  R  is  chosen  too  large. 
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Figure  13.  Optimal  probability  of  a  target  attack  versus  maximum  allowable  probability  of  a 
false  target  attack.  The  scenario  involves  a  normally-distributed  target  and  a  Poisson  held  of 
false  targets,  c  =  100,  a  =  5  [1/distance2],  a2  =  0.25  [distance2],  and  R  =  5  [distance]. 


V.  Conclusion 

Three  scenarios  were  considered  involving  an  airborne  munition  searching  for  stationary 
targets  among  a  field  of  false  targets.  Targets  were  modelled  using  uniform,  Poisson,  and 
normal  distributions.  False  targets  were  modelled  using  Poisson  distributions.  The  control 
variable  was  the  probability  of  a  target  report,  which  is  linked  to  the  probability  of  a  false 
target  report  via  the  sensor’s  Receiver  Operating  Characteristic  curve.  Optimization  prob¬ 
lems  were  formulated  that  maximize  the  probability  of  a  target  attack  while  at  the  same 
time  constraining  the  probability  of  a  false  target  attack.  Generalized  expressions  for  these 
two  probabilities  were  derived  for  all  three  scenarios.  Both  constant  and  dynamic-threshold 
solutions  were  given  for  a  sample  instance  of  each  scenario.  For  the  two  scenarios  involving 
constant  target  to  false  target  ratios,  the  optimal  dynamic  threshold  was  only  a  few  percent 
better  than  the  best  constant  threshold.  For  the  scenario  where  target  to  false  target  ratio 
changed  during  the  search,  the  optimal  dynamic  threshold  was  substantially  better  than  the 
optimal  constant  threshold. 

Future  work  will  include  parameterizing  Receiver  Operating  Characteristic  curves  by 
area  coverage  rate  and  formulating  problems  for  air  vehicles  with  multiple  warheads.  These 
generalizations  will  open  the  solution  space  to  planar  regions  and  expand  the  applications  to 
unmanned  combat  air  vehicles  and  general  sensor  craft.  As  autonomous  unmanned  vehicles 
become  more  and  more  prevalent  in  combat  environments,  more  basic  research  involving 
analytic  effectiveness  measures  (like  the  work  presented  here)  is  needed. 


23  of  30 


Appendix 


Formulae  for  probabilities  of  events  are  derived  when  parameters  vary  with  time  or  space. 
Specifically,  one  is  interested  in  the  probability  of  no  target  attacks  and  the  probability  of  no 
false  target  attacks  prior  to  a  time  t  or  radius  r.  For  example,  no  target  attacks  happen  if  no 
targets  were  encountered  or  all  targets  encountered  were  misclassihed.  One  can  denote  by 
Pc  a  general  classification  probability  which  can  assume  the  role  of  Ptr,  1  —  Ptr ,  Pftr ,  or 
1  —  Pftr  depending  on  the  application.  Events  involving  only  encounters  are  denoted  S,  and 
events  involving  classified  encounters  are  denoted  C.  The  number  of  events  occurring  in  a 
given  area  is  subscripted  along  with  the  applicable  area.  For  example,  1  classified  encounter 
in  AAi  is  denoted  C^aAi- 


A.  Uniform  Distribution 

Consider  the  case  of  searching  for  one  uniformly-distributed  target.  The  air  vehicle  flies  a 
straight  path  as  illustrated  in  Figure  1(a);  however,  its  classification  probabilities,  search 
width,  and  velocity  are  allowed  to  vary  with  time.  The  probability  of  no  encounters  (and 
hence  no  classification)  during  [0,£]  is 

P  (So, a)  =  1  -  -7-  [  w(t)v  (t)  dr  (70) 

As  JO 

where 

As  —  I  w  (r)  v  (r)  dr.  (71) 

Jo 

To  calculate  the  probability  of  exactly  one  classified  encounter  occurring  during  [0,£],  divide 
the  time  interval  into  n  short  time  periods  of  length  Ti, ...  ,Tn  such  that 

n 

!>  =  «  (72) 

i= 1 

Let  the  mean  probability  of  classification,  search  width,  and  velocity  in  the  ith  interval  be 
Pc.,  Wi ,  and  Vi  respectively.  The  incremental  area  is  given  by 

AAi  =  WiViTi  (73) 


and  the  search  area  is  given 


n 


As  =  ^2 

i 


(74) 


24  of  30 


The  probability  of  exactly  one  classified  encounter  occurring  during  [0,  t]  requires  considering 
all  the  permutations  where  one  classified  encounter  could  occur.  That  is 

P  (Ci ,a)  =  P  (Ci.AAi)  +  P  (C1Aa2)  +  •  •  •  +  P  {CiAAn)  (75) 


so 


D  tn  \  wlvlTl  D  I  D  ,  ,  WnVvTn 

1  \yi,A)  ~  - 1 - Uci  H - ; - q.,  +  •  •  •  + 


A. 


A, 


Pc"  =  4 WiViTiPCi '  (76) 


4=1 


Taking  the  limit  as  t*  — »  0,  i  =  1, . . . ,  n,  n  — >  oo  such  that  Xu=i  r*  =  ^  giyes 


T3  (Ci,A)  =  - -  Pc(t)w  (r)  v  (t)  dr. 


A 


(77) 


s  JO 


Now  calculate  the  joint  probability  of  no  encounters  and  one  classified  encounter  during  [0,  t]. 
The  two  events  are  mutually  exclusive,  thus  one  must  sum  the  two  individual  probabilities 


P  (Co, A  n  Cl, a)  =  1  -  -J-  I  W  (r)  v  (t)  dr  +  —  /  Pc  (r)  w  (r)  v  (r)  dr 


S  JO 
1  '■* 


A 


S  JO 


=  1  -  -J-  I  [1  -  Pc  (t)]  w  (r)  v  (r)  dr. 


(78) 


^  JO 


For  the  probability  of  no  target  attack  prior  to  time  t,  Eq.  (78)  is  used  with  Pc  (r)  = 
1  -  Ptr  (t)  giving 

1  /'* 

(79) 


1  " 
^  Jo 


pta  (0  =  1  -  /  pTR  (t)  w  (r)  v  (r)  dr 


B.  Poisson  Distribution 

Let  jj  (r)  be  the  rate  of  encounters  at  time  r.  Divide  the  time  interval  [0,  t]  into  n  short  time 
periods  of  length  Ti, ...  ,Tn  such  that 


Ti  =  t  (80) 

i= 1 

Let  the  mean  rate  of  occurrence  and  mean  probability  of  classification  in  the  ith  interval  be  /x, 
and  PCi  respectively.  The  standard  Poisson  probability  law  requires,  “If  an  area  is  subdivided 
into  n  subareas  and  for  i  —  1, . . . ,  n,  denotes  the  event  that  at  least  one  or  more  encounters 
occur  in  the  ith  subarea,  then,  for  any  integer  n,  Ei, . . . ,  En  are  independent  events.”3  Hence, 
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the  probability  that  exactly  jt  encounters  occur  in  the  interval  t*,  i  —  1, . . . ,  n  is 


p(h . a.)  = 


\3i 


i=  1 


En 

_  ^  i= 1  A HTi 


n 

i=  1 


(j-liTi 


30 


(81) 


The  probability  of  no  encounters  occurring  (and  hence  no  classification)  during  [0,  t]  is 


P  (Co,  a)  =  e 


—  P  Si= 1  WT> 


(82) 


Let  j  —  1-  The  probability  of  one  classified  encounter  in  the  ith  interval  is 


P(C 


Pr ... 


(83) 


The  probability  of  exactly  one  classified  encounter  occurring  during  [0,  t]  requires  considering 
all  the  permutations  where  one  classified  encounter  could  occur.  That  is 


P  (Cl ,a)  —  P  (Ci.AAl)  +  P  (Cl,AA2)  +  '  •  '  +  P  (Cl, A An) 


(84) 


So, 


P(Ci,a)  =  e 


-  p-r.'L,  P  I  e-yLwTiWV  i  ...  i  e~y:Li  wri  (^n)1  p 

^  1 !  C2  r  TC  1!  1  Cn 


=  e  sr=i  w^HiTiPci 

i= 1 


(85) 


The  probability  of  exactly  two  classified  encounters  occurring  during  [0 ,t]  requires  consider¬ 
ing  all  the  permutations  where  two  classified  encounters  could  occur.  That  is 


P(C2,a) 


P  (CgAAi  n  Ci,AA2)  +  P  (CgAAl  n  Ci;AA3)  +  •  •  •  +  P  (Cl,AAi  n  Ci,AAn)  + 
c*  (Ci)AA2  hi  C\,AAz)  +  P  (Ci,AA2  hi  CiiAA4)  +  '  '  •  +  -P  (CiiAA2  hi  Ci;AA„)  + 


+-P  (Ci,AA„_i  hi  Ci,AAn)  + 

-P  (C2,AAi)  +  P  (C2,AA2)  +  •  •  •  +  P  (C2,AA„)  (86) 
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So, 


P(C2,a)  =  e-sr=iwri(^ TjLpci 

(>2  r2)X 


(/i2r2)  p  {^nTn)  p 

*C2  1***1  -TCri 


En 

i=l  /^i7* 


P 

1!  C2 


1! 

(pPa)1 

1! 


P  (f1»7n)  p 

-rC3  1“  “T  ^  |  1  Cn 


+ 


+ 


3-  E?=1  wn  l7~n— l)  p  (^»T»)  p  , 

2  |  ^Cn_l  2  I  ^cn  ' 

(h 1T1 )  n2  |  i  if^nTn)  p2 

2!  ^Cl  '  '  "  "  2! 


g  E"=l  AtiTi 

1 


_  —  g  E,_iftr*  [2//iTiPCi  (h2T2PC2  +  •  •  •  +  UnTnPcn)  + 

2/i2 T2FC2  (/73T3Pc3  H - f-  HnTnPcn)  + 

“I”  2/^n—  lT^i— 1-fcn—j l^n^nPcn "b 

(hiri PCl)2  +  (id2T2PC2)2  H - h  (^„rnFc„)2] 

=  e-  E?=i  (E"=i  ViTiP*)2 


(87) 


In  general,  the  probability  of  exactly  j  classified  encounters  occurring  during  [0,  £]  is 

p  {Cj  a)  =  e~  sr=i  wn  (Sr=i/^P*)J  (88) 

j! 

Taking  the  limit  as  r*  — >  0,  i  =  1, . . . ,  n,  n  — >  oo  such  that  )T)”=1  t,:  =  £  yields 


P(Cj-a)  =  e~tiPT)dT 


fo  Pc  (r)  h  (t)  dr 


(89) 


Equation  (89)  can  be  used  to  determine  the  probability  of  either  no  false  target  attacks  or 
no  target  attacks  occurring  up  through  time  £.  In  the  case  of  no  false  target  attacks,  the 
possible  mutually  exclusive  events  are 

no  encounters 

1  encounter  correctly  classified 

2  encounters  correctly  classified 


oo  encounters  correctly  classified 
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Thus, 


P  (£o ,A  n  ChA  n  .  .  .  n  Coo,a)  —  P  (£o,a)  +  P  (Cl ,a)  +  •  •  •  +  P  (Coo, a) 


—  g-  fo  Ppdr  _j_  e-  f*  n(T)dT 


fo  pc  (r)  P  (r)  dr 

T\ 


+ 


_ J_  g-  fo  PT)d' 

—  g-  fo  M T)dr 


Jo  (r)  /x  (r)  dr 


CX)! 


1  ,  (fo  PdT)PT)dr)  (/o  Pc{r)fi(T)dT 

I  1  I  1***1 


=  fo  v(T)dTpfo  Pc(r)/j,(T)dT 


—  p-/o[1-Pc(T)]#1(T)dT 


(90) 


To  calculate  the  probability  of  no  false  target  attacks,  let  Pc  (r)  =  Pftr  (t)  and  //  (r)  = 
a  (r)  to  (r)  v  (r)  in  Eq.  (90)  giving 


_  g-/o[1--pFTii('r)]a('r)«)(r)?j(r)(iT 


(91) 


where  a  (r)  is  the  false  target  density.  To  calculate  the  probability  of  no  target  attacks,  let 
Pc  (r)  =  1  -  PTR  (r)  and  \i  (r)  =  P  (t)  w  (t)  v  (r)  in  Eq.  (90)  giving 

P^  (t)  =  e“  /o  [pT«(r)]/3(r)«,(r)i;(r)dr  (92) 


where  (3  (r)  is  the  target  density.  The  integrals  in  Eq.  (91)  and  Eq.  (92)  serve  as  Poisson 
parameters,  which  are  assumed  constant  in  the  literature.  These  integrals  are  “dynamic  Pois¬ 
son  parameters” .  When  formulating  optimal  control  problems,  states  are  used  to  represent 
the  integrals. 


C.  Normal  Distribution 


Consider  the  case  of  searching  for  one  normally-distributed  target  in  a  circular  disc  of  radius 
r.  Normally-distributed  refers  to  circular  normal  distribution  with  standard  deviation  a. 
The  classification  probabilities  vary  with  radius.  The  probability  of  no  encounters  (and 
hence  no  classification)  searching  a  disc  of  radius  r  from  the  origin  outward  using  concentric 
annuli  of  thickness  dp  is 


P(So,a) 


[2n  nr  1  p2  P  1  _  P2 

1-  /  /  - — 2e  2°spdpd6  =  l-  j  —e  ^ pdp. 

Jo  Jo  2707-  J0  a 


(93) 
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To  calculate  the  probability  of  one  classified  encounter,  divide  the  disc  of  area  A  into  n 
concentric  annuli  of  area  Ai, . . . ,  An  such  that 


E  -4<  =  A  (94> 

i— 1 


Each  Ai  is  calculated  using 


Ai  =  27 rpiApi 


(95) 


where  pi  is  the  radius  of  the  annulus  and  A  pi  is  the  thickness.  The  search  begins  at  the 
origin  of  the  disc  and  progresses  outward  such  that  pi  <  pt+i,  i  =  1, . . .  ,n  —  1.  Let  the 
mean  probability  of  classification  in  the  ith  annulus  be  PCi.  The  probability  of  one  classified 
encounter  during  the  sweep  of  A  requires  considering  all  the  permutations  where  one  classified 
encounter  could  occur.  That  is 


P  (Cl  ,a)  —  P  (Ci,AAi)  +  P  (Cl, A  A2)  +  •  •  •  +  P  (Cl, A  An) 


(96) 


SO 


1  _  pi  l  _  pI 

P  (C itA)  =  27rp1Ap-i— ^ PCl  H - h  2npnApn^—^e  ^ 


1  p ? 

Pc„  =  —  P-ApreA  PCt. 


(97) 


Taking  the  limit  as  A pt  — >  0,  i  =  1, . . . ,  n,  n  — >  oo  such  that  YH=\  A  =  A  yields 


P  (Ci, a)  = 


A(p)  -4  A 
— —  e  2^2  pap. 

a2 


(98) 


Now  calculate  the  joint  probability  of  no  encounters  and  one  classified  encounter  during  the 
sweep  of  A.  The  two  events  are  mutually  exclusive,  thus  one  must  sum  the  two  individual 
probabilities 


P  (£0ja  H  CljA)  =  1-1  \e  Apd,p  + 


=  1  - 


° 

[1-Pc{p)\  -A  , 

1 5 - -e  ^ pclp. 

a2 


PM  -A  i 

— —  e  2^2  pdp 
a2 


(99) 


For  the  probability  of  no  target  attack  prior  to  radius  r,  Eq.  (99)  is  used  with  Pc  (p)  = 
1  -  Ptr  (p)  giving 

rr  id  I  2 

(100) 


Pta  (r)  —  1 


J  0 


Ptr  ip)  -A  a 
- - — e  ^  pdp 


a* 
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